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The results of a preliminary partially quenched (Nf = 2) study of the light hadron spectrum, renormalization 
constants and light quark masses are presented. Numerical simulations are carried out with the LL-SSOR 
preconditioned Hybrid Monte Carlo with two degenerate dynamical fermions, using the plaquette gauge action 
and the Wilson quark action at /3 = 5.8. Finite volume effects have been investigated employing two lattice 
volumes: 16 3 x 48 and 24 3 x 48. Configurations have been generated at four values of the sea quark mass 
corresponding to Mps/Mv 0.6 0.8. 



Lattice QCD calculations of the light hadron 
spectrum and quark masses have significantly 
improved in recent years. In the quenched ap- 
proximation, i.e. ignoring quark vacuum polar- 
ization effects, the accuracy has become smaller 
than 10% and the quenching error remains then 
the main source of uncertainty. In the present 
work we explore sea quark effects in light hadron 
and quark masses, with two dynamical flavors, 
using the plaquette gauge action and the Wilson 
quark action at (3 = 5.8. 

Configurations have been generated by us- 
ing the LL-SSOR preconditioned Hybrid Monte 
Carlo |1I2| . We have implemented the Leap- 
Frog integration scheme, with trajectory length 
equal to one and time step 8t = 4 • 10~ 3 , and 
the BiCGStab inversion algorithm with iterated 
residual r = 10 -15 . The acceptance probability 
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is found to be larger than 80% and the inversion 
rounding error results to be well under control, 
being of the order of 10~ 8 . We have verified 
that reversibility is satisfied at the relative level 

of nr 11 . 

By looking at the autocorrelation times of the 
plaquette and of the pseudoscalar and vector 
two point correlation functions, we made the 
conservative choice of selecting configurations 
separated by steps of 45 trajectories. 

We investigated finite size effects by simulat- 
ing two lattice volumes: 16 3 x 48 (L ~ 1.1/m) 
and 24 3 x 48 (L ~ 1.6/m). At the smaller vol- 
ume 100 configurations have been generated at 
each of the four values of the sea hopping pa- 
rameter (k sea = 0.1535,0.1538,0.1540,0.1541), 
corresponding to Mps/My — 0.6 ~ 0.8. At the 
larger volume, where the simulation is still in 
progress, we have generated 50 configurations 
at k sea = 0.1535,0.1540 and 25 configurations 
at k sea = 0.1538,0.1541. 
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Figure 1. Degenerate (k v i = k v2 = k 3ea = 0.1540) 
pseudoscalar and vector effective masses, at the two 
simulated lattice volumes: 16 3 x 48, 24 3 x 48. 



1. Pseudoscalar and vector meson masses 

In order to study the light hadron spectrum 
and quark masses we have calculated two points 
pseudoscalar and vector correlation functions 
with both degenerate and non-degenerate va- 
lence quarks and with valence quark masses 
equal or different to the sea quark mass. 

We have tried to implement Jacobi smear- 
ing 0, but we don't find a significant im- 
provement in the final determination of hadron 
masses; the main effect is that the ground state 
can be isolated at smaller times. 

The quality of the plateaux of both pseu- 
doscalar and vector meson effective masses is 
shown in FigQ] One can also see from the 
plot that finite volume effects are found to be 
at the level of ~ 20% on the smaller volume 
V = 1 6 3 x 48. Such large effects are perhaps un- 
expected since they should scale as exp (—ML), 
where ML ~ 4.4 for pseudoscalars and ~ 5.8 
for vectors in our simulation on the small lat- 
tice. We believe that this point requires further 
investigations. For the time being, the results 
we present below follow from the analysis per- 
formed at the larger volume V — 24 3 x 48. We 
observe a smooth enhancement of finite volume 
effects as the sea quark mass decreases. In Fig|5] 
left (right) we show the pseudoscalar(vector) 
meson masses as a function of the valence quark 
mass (l/k v — l/2(l/fc„i + l/fc„ 2 ))- Different 
symbols refer to different k sea values. Within 
the statistical accuracy, meson masses are found 
to be linear in the valence quark mass, while 
the dependence on the sea quark mass is less 
clear. In particular both pseudoscalar and vec- 
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Figure 2. Pseudoscalar squared(left) and vec- 
tor(right) meson masses as a function of the valence 
quark mass. Different symbols refer to different val- 
ues of the sea quark mass. 

tor meson masses obtained at k sea = 0.1538 are 
smaller than those at k sea = 0.1540, although 
compatible within the errors. As a consequence, 
the dependence on the sea quark mass is not 
easily taken into account in the chiral extrapo- 
lations. 

In order to determine the lattice spacing and 
the values of the bare light quark masses we 
fit the vector meson masses linearly in the 
squared pseudoscalar masses. Following the 
method of "lattice physical planes" 0], we set 
the ratio Mk/Mk* equal to its experimental 
value, thus obtaining for the inverse lattice spac- 
ing the estimate a -1 = 3.0(I)GeV. For the 
pseudoscalar and vector meson masses we find 
M„ = I43(3)MeV, M p = 805(19)MeV (with 
M„/M p fixed to its experimental value) and 
M$ = 983(1 9)MeV. 

2. Renormalization Constants 

The determination of the renormaliza- 
tion constants has been performed non- 
perturbatively by using the RI-MOM method*. 
We have considered the bilinear quark opera- 
tors O r = qTq with T — S,P, V, A, T standing 
respectively for /, 75, 7 P , 7^75, and a^. We also 
present the result for the renormalization con- 
stant of the quark field, Z q . 

Finite volume effects are found to be smaller 
than 5% for all the renormalization constants 
and all the values of the sea quark mass. 

The dependence of the renormalization con- 
stants on the sea quark mass is weak, as one 
can see form Fig|3 where the scale independent 
combination Zp GI |Sj is shown as a function of 
the renormalization scale. 



'For a more detailed explanation see ref. |5] and refer- 
ences therein. 
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Figure 3. Zp at different values of the sea quark 
mass, as a function of the renormalization scale 
(afi) 2 . At large scales discretization effects intro- 
duce a smooth linear dependence on (an) 2 . 

The results for the renormalization constants, 
extrapolated to the chiral limit, in the RI-MOM 
scheme at the scale \i = a -1 , read 

Z s = 0.72(1) , Z P = 0.55(1) , Z v = 0.69(2) , 
Z A = 0.79(2) , Z T = 0.77(5) , Z q = 0.79(1) . (1) 



3. Light quark masses 

In order to determine the values of quark 
masses we used, as in ref. |H|, the standard 
procedures based on the vector (V) and axial- 
vector (A) Ward Identities(WI). 

The VWI connects the quark mass renormal- 
ization constant to the scalar density one. The 
corresponding mass definition is given by 



(2) 



where k q is the Wilson hopping parameter 
and k cr is its critical value, corresponding to 
Mp s (k v — k sea = k cr ) — 0. The AWI, instead, 
relates the quark mass renormalization constant 
to the renormalization constant of the axial and 
pseudoscalar operators, leading to the following 
expression for the renormalized quark mass 
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(3) 



In order to get the physical values of light 
quark masses, we study the dependence of the 
squared pseudoscalar masses on simulated va- 
lence and sea quark masses. The observed be- 
haviour on the sea quark mass doesn't allow 



us to extract from the fit the dependence on 
(am sea ), and the dependence on the valence 
quark mass is found to be linear. We perform, 
therefore, a fit to the form 
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and a similar expression for the VWI quark 
mass. The constant term A in eq. |@J is due 
to 0(a)— discretization effects and it is only 
present in the AWI case. For the VWI case 
these effects are automatically included in the 
determination of the critical hopping parame- 
ter k cr , implying A = 0. The physical values 
of the average up/down (mi) and of the strange 
(m s ) quark masses are then obtained by substi- 
tuting the experimental pion and kaon masses 
on the l.h.s. of eq.0]and the value of the lattice 
spacing (a' 1 = 3.0(l)GeV). 

Quark mass values are converted from RI- 
MOM at the renormalization scale \i = l/a to 
MS at fi = 2GeV by using RG improved per- 
turbation theory at 4— loop accuracy j7 • The 
preliminary results read 



ywi 



AWI 



4.8(5)MeV , m 



VWI 



lll(6)MeV, 



= 4.5(5)MeV,mf wl = 103(9)MeV. (5) 

Simulations at other values of /3, k sea and larger 
volumes are required, in order to study the con- 
tinuum limit, the dependence on the sea quark 
mass and finite volume effects. 

It is a pleasure to thank the Lattice-2004 or- 
ganizers for the pleasant and stimulating con- 
ference realized at the Fermilab. 
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